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The  coolant  heat  sinks  in  thermoelectric  generators  (TEG)  play  an  important  role  in  order  to  power  genera¬ 
tion  in  the  energy  systems.  This  paper  explores  the  effective  pumping  power  required  for  the  TEGs  cooling  at 
five  temperature  difference  of  the  hot  and  cold  sides  of  the  TEG.  In  addition,  the  temperature  distribution  and 
the  pressure  drop  in  sample  microchannels  are  considered  at  four  sample  coolant  flow  rates.  The  heat  sink 
contains  twenty  plate-fin  microchannels  with  hydraulic  diameter  equal  to  0.93  mm.  The  experimental 
results  show  that  there  is  a  unique  flow  rate  that  gives  maximum  net-power  in  the  system  at  the  each  tem¬ 
perature  difference. 
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1.  Introduction 

In  the  power  generation  systems,  a  key  factor  is  the  optimization  of 
the  systems  design,  together  with  its  heat  source  and  heat  sink.  In  the 
case  of  the  thermoelectric  generators  (TEG),  for  increasing  the  convec¬ 
tive  surface  area,  in  order  to  provide  high  density  of  heat  dissipation  at 
the  surface  of  the  structure,  microscale  heat  transfer  systems  can  en¬ 
hance  the  thermal  coupling  to  the  hot  and  cold  reservoirs  [  1  ].  Therefore, 
the  challenge  is  to  design  an  effective  heat  exchanger  within  microelec¬ 
tronic  dimension  restrictions  [2].  Microscale  single-phase  heat  transfer 
has  been  widely  used  in  industrial  and  scientific  applications  [3]. 
Using  microchannel  heat  sinks  provides  low  weight  and  compact  ener¬ 
gy  systems,  compared  to  the  traditional  macroscale  heat  sinks,  and  in¬ 
creases  modularity.  In  contrast  to  macrochannels,  a  reduced  flow  rate 
in  the  microchannel  heat  sink  is  sufficient  to  maintain  the  same  average 
temperature  difference  between  the  hot  and  cold  sides  of  the  TEGs. 
Using  microchannel  heat  sinks  also  increases  the  Nusselt  number  in 
the  channels  [4]. 

In  microchannels,  the  convective  heat  transfer  increases  at  the 
high  relative  roughness  of  the  channel  walls,  so  that  the  critical  Reyn¬ 
olds  number  changes  with  the  wall  roughness  in  a  different  way  com¬ 
pared  to  the  typical  macrochannels.  There  are  observations  of  earlier 
transition  from  laminar  to  turbulent  flow  regime  in  microchannels 
compared  to  the  macrochannels  theory  [5,6].  The  earlier  transition 
happens  at  lower  Reynolds  number  when  the  hydraulic  diameter  of 
the  microchannel  decreases  [7,8].  The  studies  indicate  that,  the  geo¬ 
metric  configuration  of  the  microchannel  heat  sinks  has  a  critical  ef¬ 
fect  on  the  convective  heat  transfer  of  the  laminar  flow  in  the  heat 
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sink  [9].  In  order  to  achieve  overall  heat  transfer  enhancement  a  rectan¬ 
gular  microchannel  is  the  best  shape,  and  its  heat  transfer  coefficient  is 
the  highest  amongst  trapezoidal  and  triangular  shaped  microchannels 
[10]. 

A  comparison  of  thermal  efficiencies  is  conducted  with  triangular, 
rectangular  and  trapezoidal  microchannels  with  the  same  hydraulic  di¬ 
ameter  by  Chen  et  al.  [11].  Their  results  show  that  the  cross-sectional 
shape  has  a  significant  effect  on  the  temperature  distribution  of  heat 
sinks,  and  the  lowest  pumping  power  is  required  for  the  triangular 
microchannel  heat  sink.  Kroeker  et  al.  [12]  found  that,  on  the  basis  of 
equal  hydraulic  diameter  and  equal  Reynolds  number,  a  rectangular 
channel  has  less  thermal  resistance  compared  with  heat  sinks  with  cir¬ 
cular  channel.  Additionally,  the  laminar  friction  factor  or  flow  resistance 
reaches  a  minimum  value  as  the  channel  aspect  ratio  approaches  0.5  [9]. 
As  reported  by  Tuckerman  and  Pease  [13],  when  the  channel  width  and 
the  fin  thickness  are  equal,  the  heat  sink  gives  the  minimum  thermal  re¬ 
sistance.  The  inlet/outlet  plenum  effect  gives  a  non-uniform  velocity 
and  flow  rate  distribution  in  each  channel  under  a  given  pressure 
drop  in  the  heat  sink,  so  that  a  non-uniform  temperature  distribution 
happens  in  the  heat  sink  [14].  The  heat  sinks  with  the  vertical  coolant 
supply  and  collection  via  inlet  and  outlet  ports  of  the  heat  sink  provide 
better  uniformities  in  temperature  and  velocity,  so  that  this  type  of  heat 
sink  can  make  better  performance  due  to  smaller  thermal  resistance 
among  the  studied  heat  sinks. 

The  heat  exchanger  designs  have  mostly  disconnected  to  the  higher 
performance  TEGs.  The  development  work  mostly  focused  on  thermo¬ 
electric  materials  required  a  significant  amount  of  engineering  para¬ 
metric  analysis.  Since,  the  ZT  of  the  thermoelectric  materials  is  not  the 
only  factor  to  improve  the  output  power  of  the  system,  the  whole 
parts  of  the  energy  conversion  system,  such  as  thermal  contacts  of  the 
hot  and  cold  heat  exchangers,  needed  to  be  considered.  The  thermal 
resistances  of  the  heat  exchangers  have  a  strong  influence  on  the 
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Nomenclature 

A  TEG  leg  area,  m2 

Cp  specific  heat  of  water,  4178  J/kg.K 

D  diameter,  m 

Dh  hydraulic  diameter  of  the  channel,  m 

H  height,  m 

I  current,  A 

L  length,  m 

m  mass  flow  rate,  kg/s 

P  power,  W 

p  pressure,  Pa 

Pr  Prandtl  number 

Re  average  Reynolds  number 

Q.  heat  absorbed  by  heat  sink,  W 

T  temperature,  K 

t  heat  sink  wall  thickness,  m 

u  coolant  flow  velocity  in  the  microchannel,  m/s 

V  voltage,  V 

V  volumetric  flow  rate,  m3/s 

w  width,  m 

xh  thermal  entry  length,  m 


Greek  symbols 

a  Seebeck  coefficients,  V/K 

A  value  difference 

/d  water  dynamic  viscosity,  N.s/m2 

p  water  density,  kg/m3 


Subscripts 

av  average 

ch  microchannel 

hs  heat  sink 

i  inlet 

n  net 

o  outlet 

p  pump 

pi  plenum 

teg  thermoelectric  generator 


electric  power  produced  by  the  TEG.  Martinez  et  al.  [15]  optimized  the 
heat  exchangers  of  a  TEG  in  order  to  maximize  the  electric  power  gen¬ 
erated  by  the  TEG.  In  addition,  in  order  to  decrease  the  coolant  pumping 
power  in  the  TEG  systems,  an  effective  design  of  the  microchannel  heat 
sinks  is  proposed  and  implemented  in  a  three-dimensional  TEG  model 
[16].  The  small  thermal  conductivity  of  the  thermoelectric  materials 
causes  the  temperature  difference  between  the  cold  and  the  hot  sur¬ 
faces  of  TEG  not  to  change  remarkably  when  the  coolant  flow  rate  varies 
in  the  heat  sink.  Due  to  the  Seebeck  effect  in  the  TEG,  which  transfers  a 
certain  amount  of  energy  from  high-temperature  fluid  to  electricity,  the 
linear  variations  in  temperature  of  the  fluids  in  the  TEG  system  are  dif¬ 
ferent  from  the  logarithmic  variation  case  in  the  ordinary  parallel-plate 
heat  exchanger  [17].  Therefore,  the  heat  flow  from  the  hot  source  is  not 
equal  to  the  heat  flow  to  the  low  temperature  fluid.  Further  analysis  of 
the  influence  of  coolant  flow  rate,  heat  exchanger  geometry,  and  fluid 
inlet  temperature  on  the  power  generation  by  the  TEGs  is  done  by 
Esartea  et  al.  [18]. 

Nonetheless,  when  it  comes  to  real-world  design  of  thermoelectric 
systems  for  direct  thermal  to  electricity  conversion,  a  narrow  discussion 
of  the  TEGs  integrated  to  heat  exchangers  in  micro  scale  is  still  lacking. 
Wojtas  et  al.  [19]  reported  a  complete  integration  of  the  microfluidic 
heat  transfer  system  with  a  micro-TEG,  where  the  thermoelectric  Seebeck 


voltage  is  measured  as  a  function  of  the  temperature  gradient  and  applied 
fluid  flow  rate.  Their  demonstration  enables  us  to  increase  the  output 
power  performance  of  the  micro-TEGs.  The  thermoelectric  Seebeck  volt¬ 
age  is  measured  as  a  function  of  the  applied  flow  rate  and  temperature 
gradient.  In  this  work,  a  micro  plate-fin  heat  sink  that  is  made  of  alumi¬ 
num  is  applied  to  a  TEG.  The  purpose  is  considering  the  power  generation 
versus  the  pumping  power,  which  is  produced  by  the  pressure  loss  of  the 
coolant  fluid  in  the  system.  The  particular  focus  of  this  experimental 
study  is,  exploring  the  optimum  coolant  flow  rate  that  results  the  maxi¬ 
mum  net-power  in  the  system,  at  a  wide  average  temperature  difference 
range  of  the  hot  and  cold  sides  of  the  TEG  ( ATteg,av)-  By  considering  the  av¬ 
erage  temperature  of  the  hot  surfaces  of  the  TEG  legs,  these  optimum 
points  are  suggested  as  the  optimum  applied  pressure  drop  in  the 
microchannel  heat  sink  by  Rezania  and  Rosendahl  [20].  In  addition,  the 
temperature  distribution  and  pressure  loss  in  some  microchannels  with 
variation  of  the  coolant  flow  rate  and  heat  flux  are  reported.  Fig.  1  illus¬ 
trates  the  schematic,  thermocouple  and  pressure  gages  location,  and  the 
setup  of  the  experiments. 

2.  Theoretical  background 

The  heat  transfer  rate  removed  by  the  coolant  flow  in  the  heat  sink 
is  [21]: 

Q  =  mCp(T0-Ti).  (1) 

The  pumping  power  to  circulate  the  coolant  flow  in  the  heat  sink 
is  related  to  the  pressure  drop  and  the  volumetric  flow  rate,  and  can 
be  calculated  as  follows: 

Ppump  =  VAp  =  V(Po-Pi).  (2) 

The  thermal  entry  length  of  the  internal  flow  at  the  laminar  re¬ 
gime  is  approximated  by  [22]: 


xh^0.05DhRePr. 


(3) 


where  Dh  is  the  hydraulic  diameter  of  a  rectangular  cross  section 
shape  channel  defined  as  follows: 


Dh 


2  wH 
w  +  H' 


(4) 


The  Reynolds  number  in  the  channels  is  [23]: 
Re  =  ^. 


(5) 


In  this  study,  the  effect  of  the  temperature  variation  on  the 
thermo-fluid  properties  of  the  coolant  flow  is  taken  into  account.  In 
the  TEGs,  if  the  hot  and  cold  junctions  are  maintained  at  different 
temperatures,  an  open-circuit  electromotive  force  develops  between 
them,  and  is  given  by  V=  oc  AT,  which  defines  the  differential  Seebeck 
coefficient  between  the  two  sides  of  the  elements.  Therefore,  the 
power  generation  is  as  follows  [24]: 

Pteg  =  IV  =  oc2a  X  ATteg2  Xp;  (6) 

Lteg 


where  oc  2a  is  the  power  factor  of  the  thermoelectric  materials.  We 
define  the  net-power  in  the  system  as  follows: 

P  —  P  — P 

1  n  1  teg  1  p- 


(6) 
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Fig.  1.  The  geometrical  configuration  of  the  system,  a.  System  schematic,  b.  Thermocouples  location,  c.  Pressure  gages  location,  d.  Setup  of  the  experiments. 


3.  Experimental  procedure 

The  fabricated  microchannel  heat  sink  for  this  work  contains  twenty 
rectangular  cross  section  shape  microchannels.  The  inlet  and  outlet  ple¬ 
nums  of  the  heat  sink  open  on  the  top  cover  plate,  and  the  coolant  flow 
passes  through  a  U-shape  path  in  the  heat  sink.  Therefore,  it  is  supplied 
and  leaves  the  heat  sink  vertically.  According  to  [14],  this  type  of  the 
plenum  makes  a  better  uniform  velocity  distribution  in  the  heat  sink, 
so  that  it  provides  better  surface  temperature  uniformity  in  the  heat 
sink.  Base  on  this  uniform  mass  flow  rate  distribution  in  the  channels, 
an  average  Reynolds  number  can  be  valid  for  the  channel  through  entire 
the  heat  sink.  The  imposed  coolant  flow  rates  in  this  study  generate 
Reynolds  number  between  Re  =  63  and  Re  =  1473.  This  range  of  the 
Reynolds  number  is  lower  than  the  reported  critical  number  by  Li  et 
al.  [25].  They  found  that  the  transition  to  the  turbulence  regime  began 
near  Re  =  1 535  in  microchannels,  which  is  lower  than  critical  Reynolds 
number  predicted  by  the  classical  theory.  Therefore,  in  this  study,  the 
maximum  imposed  flow  rate  in  the  heat  sink  is  1.51/min.  The  estimated 
thermal  entry  length  varies  between  Xh  =  2.06  mm  for  the  minimum 
average  Reynolds  number  (Re  =  63)  and  xh  =  48  mm  for  the  maxi¬ 
mum  average  Reynolds  number  (Re  =  1473).  Although,  the  thermal 
entry  length  is  short  at  low  flow  rate  but  as  the  flow  rate  increases, 
the  effect  of  the  thermal  entry  length  cannot  be  negligible. 

In  order  to  measure  the  temperature  variation  of  the  heat  sink  and 
the  temperature  difference  of  the  hot  and  cold  sides  of  the  TEG,  T-type 
thermocouples  are  used.  The  thermocouples  for  measuring  the  temper¬ 
ature  of  the  coolant  fluid  in  the  microchannels  are  installed  in  the  solid 
region  of  the  heat  sink  substrate  plate,  but  close  to  the  fluid  region.  Since 
the  height  of  the  microchannels  is  small  compared  to  the  channel 
width,  the  thermocouples  represent  an  accurate  temperature  of  the 


bulk  fluid  in  the  microchannels.  By  using  eight  pressure  gages  at  the 
inlet  and  outlet  of  the  heat  sink  and  microchannels  #2,  #6,  and  #10 
the  pressure  drop  in  the  heat  sink  and  microchannels  are  measured. 
Fig.  lb  and  c  illustrates  the  thermocouple  and  pressure  gage  locations 
in  the  heat  sink,  respectively.  The  model  of  the  TEG  used  for  this  exper¬ 
iment  with  dimension  of  56  x  56  mm  is  G2-56-0375  made  by  Tellurex. 
The  voltage  and  current  of  the  TEG  are  regulated  to  give  maximum 
power  at  any  ATteg. 

The  data  of  this  work  are  recorded  at  thermally  steady  state  in  the 
system  for  five  Artegiav  (10,20,40,60,  and  80  K)  and  four  flow  rates 
(0.07,0.5, 1.0,  and  1.5  1/min)  in  the  heat  sink.  The  ultrasonic  flow  sen¬ 
sor  for  water  continuous  measurement  is  from  Burkert  with  accuracy 
of  0.01  %.  The  inlet  coolant  flow  temperature  is  fixed  at  301.5  K  and 
the  ambient  temperature  is  approximately  constant  at  296  K  through 
the  experiments.  Water  is  used  as  the  coolant  fluid  in  the  heat  sink. 
An  appropriate  Labview  program  is  applied  to  record  the  data  of  the 
measurement.  The  geometry  details  of  the  microchannel  heat  sink  is 
mentioned  in  Table  1. 


4.  Results 

As  Fig.  2  reveals,  the  pumping  power  increases  with  the  coolant 
flow  rate  with  sharper  slope  compared  to  the  increasing  of  the  output 
power  of  the  TEG.  Although,  the  imposed  heat  flux  on  the  hot  side  of 
the  TEG  is  constant  at  each  stage  of  the  test,  when  the  coolant  flow 
rate  changes,  the  output  power  increases  with  the  flow  rate.  According 
to  a  report  by  Crane  et  al.  [26],  at  a  constant  temperature  difference  of 
the  hot  and  cold  sides  of  the  TEG  (ATteg),  the  efficiency  of  the  TEG  sys¬ 
tems  integrated  with  the  heat  exchangers  increases  when  the  TEG 


Table  1 

The  geometry  details  of  the  microchannel  heat  sink. 


Parameter 

Hch(pm) 

wch  (pm) 

t  (pm) 

Dh(pm) 

tch(mm) 

whs  (mm) 

Tpi  (mm) 

f^pi,i/o  (mm) 

Dimension 

700 

1400 

1400 

933 

56 

56 

50 

5.0 
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coolant  flow  rate  (1/min) 

Fig.  2.  The  power  generation  in  the  TEG  versus  the  pumping  power. 


Fig.  3.  The  A Tteg  variation  with  the  coolant  flow  rate  at  fixed  heat  flux. 


Fig.  4.  The  variation  of  the  inlet/outlet  temperature  difference  in  microchannel  #6  with 
the  coolant  flow  rate. 

cold  side  temperature  decreases.  Therefore,  increasing  the  coolant  flow 
rate  at  constant  heat  flux  can  result  to  higher  output  power. 

In  addition,  as  the  coolant  flow  rate  in  the  heat  sink  increases,  the 
average  coolant  flow  temperature  reduces  because,  the  average  flow 
velocity  in  the  channels  raises.  Therefore,  the  heat  removal  capacity  of 
the  flow  increases,  and  the  hot  side  temperature  of  the  TEG  reduces. 
When  the  hot  side  temperature  decreases  the  heat  loss  rate  through 
the  thermal  insulation  of  the  system  reduces,  so  that  the  higher  heat 
flux  crosses  through  the  TEG,  and  therefore  the  output  power  increase. 
However,  the  plots  in  Fig.  2  are  shown  based  on  the  constant  ATtegiav, 
the  ATteg  varies  slowly  with  the  coolant  flow  rate.  This  variation  can 


Fig.  5.  The  variation  of  the  pressure  loss  in  microchannels  #2,  #6,  and  #10  with  the 
coolant  flow  rate. 

be  noticed  in  Fig.  3  for  a  sample  Artegiav.  As  is  reported  by  Wojtas  et  al. 
[19],  the  reason  of  the  small  variation  of  the  Arteg  lies  in  the  relatively 
small  impact  decreasing  of  the  thermal  resistance  of  the  heat  sink  com¬ 
pared  to  the  absolute  value. 

As  is  expected,  the  inlet/outlet  temperature  difference  of  the  cool¬ 
ant  fluid  in  the  microchannels  decreases  with  the  thermal  resistance 
of  the  heat  sink,  due  to  rising  of  the  flow  rate.  In  addition  for  sample 
microchannel  #6,  Fig.  4  shows  that  the  temperature  difference  of  the 
coolant  fluid  increases  with  the  Arteg,av.  The  inlet/outlet  arrangement 
effect  of  the  plenums  on  the  coolant  flow  distribution  is  shown  in 
Fig.  5.  The  pressure  loss  in  the  middle  channels  reveals  linear  behav¬ 
ior  when  the  flow  rate  varies.  But  the  pressure  drop  variation  in 
microchannel  #2,  which  is  fabricated  at  the  side  of  the  heat  sink, 
has  increased  exponentially  with  the  flow  rate.  Therefore,  the  dis¬ 
tributed  rate  of  the  flow  is  higher  in  the  microchannels  at  the  sides 
of  the  heat  sink  than  the  flow  in  the  microchannels  in  the  middle 
of  the  heat  sink,  as  is  shown  by  Chein  and  Chen  [14]. 

It  can  be  seen  in  Fig.  2  that  higher  ATtegiav  and  coolant  flow  rate  pro¬ 
duce  higher  output  power  in  the  TEG.  On  the  other  hand,  the  pumping 
power  increases  with  the  rising  of  the  flow  rate  in  the  heat  sink  and 
may  decrease  the  net-power  generation  from  the  system.  Therefore, 
there  should  be  a  unique  optimum  flow  rate  at  every  ATtegAV  that  gives 
maximum  net-power.  These  optimum  points  are  revealed  in  Table  2. 
The  flow  rate  equal  to  0.5  1/min  gives  the  maximum  net-power  among 
the  other  tested  flow  rates  in  this  work  for  the  ATtegiav>40  K.  At  the 
smaller  ATteg<avl  lower  flow  rate  can  produce  maximum  net-power.  It  is 
expected  that  when  the  ATteg<av  increases,  the  optimum  flow  rate  that  re¬ 
sults  maximum  net-power  shifts  to  higher  value. 

5.  Conclusions 

A  microchannel  heat  sink  is  applied  to  a  TEG  in  order  to  make  a  light 
and  compact  energy  system.  The  power  generated  versus  the  coolant 
pumping  power  is  studied  in  the  TEG.  Additionally,  the  optimum  cool¬ 
ant  flow  rates  that  result  in  the  maximum  net-power  in  the  system 
are  explored.  The  results  show  that  there  is  a  unique  coolant  flow  rate 
at  any  ATteg>av  that  makes  maximum  net-power  in  the  system.  The 
value  of  the  optimum  flow  rate  increases  when  the  ATtegiav  increases. 


Table  2 

The  net-power  generation  at  the  studied  coolant  flow  rates. 


Flow  rate  (1/min) 

ATtegav 

10  K 

20  K 

40  K 

60  K 

80  K 

0.07 

0.035 

0.164 

0.562 

1.282 

1.890 

0.5 

0.018 

0.146 

0.632 

1.320 

2.026 

1.0 

-0.103 

0.027 

0.551 

1.263 

1.996 

1.5 

-0.406 

-0.275 

0.250 

0.970 

1.741 
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Results  also  reveal  that  the  temperature  difference  of  the  hot  and  cold 
sides  of  the  TEG  increase  slowly  with  the  flow  rate  at  a  constant  im¬ 
posed  heat  flux.  The  reason  of  this  small  variation  lies  in  the  relatively 
small  impact  decreasing  of  the  thermal  resistance  of  the  heat  sink. 
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